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Abstract 



We investigate the three-dimensional 0(2) model near the critical point by Monte 
Carlo simulations and calculate the major universal amplitude ratios of the model. 
The ratio Uq = A + / A~ is determined directly from the specific heat data at zero 
magnetic field. The data do not, however, allow to extract an accurate estimate 
for a. Instead, we establish a strong correlation of Uq with the value of a used in 
the fit. This numerical a-dependence is given by A + /A~ = 1 — 4.20(5)« + 0(a 2 ). 
For the special a-values used in other calculations we find full agreement with the 
corresponding ratio values, e. g. that of the shuttle experiment with liquid helium. 
On the critical isochore we obtain the ratio = 0.293(9), and on the critical line 

the ratio Ct/Cl = 1-957(10) for the amplitudes of the transverse and longitudinal 
correlation lengths. These two ratios are independent of the used a or //-values. 
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1 Introduction 



In quantum field theory and condensed matter physics O(N) symmetric vector mod- 
els play an essential part, because they are representatives of universality classes for 
many physical systems. The universal properties of the O(N) models - the critical 
exponents and amplitude ratios, which describe the critical phenomena - are there- 
fore of considerable importance. In three dimensions the case N = 2 is a special one: 
it is the first vector model (with increasing N) showing Goldstone effects, and the 
exponent a, which controls the critical behaviour of the specific heat, is very close 
to zero. In fact, if one plots a versus N, as determined by field theory methods 
[T|]-[U, then the function is approximately linear in N and becomes negative just 
below N = 2. The proximity of a to zero made it also difficult to determine the 
type of the singularity for the specific heat in real systems. Indeed, for the lambda 
transition of helium a nearly logarithmic singularity (corresponding to a = 0) was 
first measured pj and a similar behaviour was found at the gas-liquid critical point 
0. However, with the nowadays reached experimental precision, especially that of 
the spectacular shuttle experiment with liquid helium 0, |§ there is no doubt that 
the critical exponent a is very small, but non-zero, and because it is negative the 
peak of the specific heat is finite. 

In this paper we calculate, among others quantities, the specific heat from Monte 
Carlo simulations. The determination of a from these data poses, as we shall see, 
similar problems as in experiments. Of course, there is only one value of a for the 
3d 0(2)-universality class, but it is unclear what the correct value is (see e. g. the 
survey in Table 19 of Ref. ||). We therefore pursue the strategy to calculate the 
universal ratios from our data for different a-values in the range where the actual 
value most probably is. The strongest dependence on the used a is expected for 
fits involving the universal amplitude ratio A + / A~ of the specific heat. The same is 
true for all theoretical determinations [IIJ [TTJ] of this ratio. Apart from A + j A~ we 
derive from our simulations other universal quantities and amplitude ratios, which 
characterize the 0(2)-universality class in three dimensions. 

The model which we investigate is the standard 0(2)-invariant nonlinear cr-model 
(or XY model), which is defined by 

(3H = -J ^ $x-$y ~ H- J2$ x - (1) 

<x,y> x 

Here x and y are the nearest-neighbour sites on a three-dimensional hypercubic 
lattice, <p x is a 2-component unit vector at site x and H is the external magnetic 
field. We consider the coupling constant J as inverse temperature, that is J = 1/T. 
Instead of fixing the length of the spin vectors (p x to 1 we could have introduced an 
additional term J^xl^x + -M0ir — I) 2 ] on the right hand side of the last equation. By 
choosing an appropriate A value ]12) it is then possible to eliminate leading order 
corrections to scaling. As it will turn out, these corrections are negligibile in the 
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energy density and marginal in the specific heat also with the Hamiltonian from Eq. 
([I]). Moreover, we want to combine amplitudes obtained from former simulations 
at non-zero magnetic field Jl3| using the same Hamiltonian with the amplitudes we 



determine now in order to calculate universal ratios. 

As long as H = | H | is non-zero one can decompose the spin vector X into a 
longitudinal (parallel to the magnetic field H) and a transverse component 

4 = <fe + ft , with e H = H/H . (2) 

The order parameter of the system, the magnetization M, is then the expectation 
value of the lattice average 0" of the longitudinal spin component 

M = <^X>I> = (0">. (3) 

X 

Here, V = L 3 and L is the number of lattice points per direction. There are two 
types of susceptibilities. The longitudinal susceptibility is defined as usual by the 
derivative of the magnetization, whereas the transverse susceptibility corresponds 
to the fluctuation of the lattice average (f) 1 - of the transverse spin component 

XL = H = V((^)-M 2 ), (4) 

XT = V(ft 2 >. (5) 
The total magnetic susceptibility is 

X = Xl + Xt ■ (6) 

At zero magnetic field, H = 0, there is no longer a preferred direction and the 
lattice average of the spins 

X 

will have a vanishing expectation value on all finite lattices, (<fi) = 0; the longitudinal 
and transverse susceptibilities become equal for T > T c and diverge both for T <T C 
because of the Goldstone modes [13]. Nevertheless we can use <fr to define the total 



susceptibility and the Binder cumulant by 

X = V($'), (8) 
* = (9) 

For T > T c we have \ = 2\l = ?<Xt- We approximate the order parameter M for 
H = by ID 

M ~ . (10) 
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On finite lattices the magnetization of Eq. ( JLOD approaches the infinite volume limit 
from above, whereas M as defined by Eq. (|3|) for H ^ reaches the thermodynamic 
limit from below. 

In our zero field simulations we want to measure three further observables: the 
energy density, the specific heat and the correlation length. The energy of a spin 
configuration is simply 

e = - j2 i* ■ iv > ( n ) 

<x,y> 



and the energy density e is then 



For the specific heat C we obtain 



(E)/V. (12) 



C = § = y((E 2 )-(Ef) . (13) 



= ( M 2f _ lT , ) , (14) 



The second moment correlation length is calculated from the formula 

X/F-1 
^4sin 2 (7r/L) 

where F is the Fourier transform of the correlation function at momentum = 
27re At / L, and a unit vector in one of the three directions 

F = ^(lE ex P(^)^.| 2 ) • (15) 

X 

In the simulations we compute F from an average over all three directions. Strictly 
speaking, Eq. fll4l) can only serve as a definition of the correlation length for T > T c , 
because the exponential correlation length diverges for H —>■ and T < T c . Instead 
it is possible to introduce a transverse correlation length £t on the coexistence line 
[T5| , which is connected to the so-called stiffness constant p s for d = 3 by 

£ T = p; 1 for H = 0, T <T C . (16) 

We explain later how to calculate p s . For H ^ there are two exponential correlation 
lengths, a transverse and a longitudinal one Their second moment forms 
may be computed again from Eq. fll4]) by replacing x and F with their respective 



transverse or longitudinal counterparts. 

The rest of the paper is organized as follows. First we discuss the critical be- 
haviour of the observables and define the universal amplitude ratios, which we want 
to determine. In Section 3 we describe our simulations at H — 0, the results for 
the Binder cumulant, the critical point and the correlation length. Then we analyse 
the data for the energy and the specific heat. In Section 4 we discuss as an alter- 
native the calculation of A + /A~ from the equation of state, which was obtained 
from non-zero field simulations. The following Section 5 serves to find the specific 
heat and the correlation lengths at T c , as well as the stiffness constant, from H ^ 
simulations. We close with a summary of the ratios and the conclusions. 
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2 Critical Behaviour 



In the thermodynamic limit (V — > oo) the observables show power law behaviour 
close to T c . It is described by critical amplitudes and exponents of the reduced 
temperature t = (T — T c )/T c . We note that we use here another definition of t than 
in Ref. [I3|. We will mention this point again later. The scaling laws at H = are 
for: 

the magnetization 

M = B(-tf for t < , (17) 
the longitudinal susceptibility 

Xl = C + r~< for t > , (18) 

the transverse correlation length 

6r = fort<0, (19) 

the correlation length 

f = f for i > , (20) 
for t — > ±0 the energy density 

e = e ns + T c t(c ns + —^—-\t\-A , (21) 
\ a(l — a) J 

and the specific heat 

C = C ns + — \t\~ a . (22) 
a 

The specific heat and the energy density contain non-singular terms C ns and e ns , 
which are due to derivatives of the analytic part f ns of the free energy density They 
are the values of the specific heat and energy density at T c . With our definition for 
the specific heat amplitudes we have already singled out their main a- dependencies, 
the remaining factors A ± are only moderately varying with a. 

On the critical line T = T c or t = we have for H > the scaling laws 

M = d c H 1/s or H = D C M S , (23) 

and for the longitudinal and transverse correlation lengths £l,t 

&,t = e L , T H~ Uc , v c = v/P6 . (24) 

The specific heat scales as 

C = C ns + ^H- a ° , a c = a/ (35 . (25) 
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We assume the following hyperscaling relations among the critical exponents to be 
valid 

2- a = du, 7 = (3(6 -I), dp = (3(1 + 5). (26) 

As a consequence only two critical exponents are independent . Because of the hyper- 
scaling relations and the already implicitly assumed equality of the critical exponents 
above and below T c one can construct a multitude of universal amplitude ratios Jl5| 
(see also the discussion in Ref. |9|]). The following list of ratios contains those which 
we want to determine here 

U Q = A+/A- , = e/^r , (27) 

R + = (A + y/ d e , Rj = (A'f'% , (28) 
R x = C+DcB 6 - 1 , R c = A+C+/B 2 , (29) 

and 

R A = A c D^ 1+a ^B- 2 ^ , Ql = (&/tV /v C + /d c (l/5+l) . (30) 
One of the ratios, R x , was already calculated by us from non-zero magnetic field 



simulations [JT3J , using the exponents of Ref. |T2|]. We found 

R x = 1.356(4) . (31) 

In order to normalize the equation of state, the temperature and the magnetic field 
in the same paper, we had computed the critical amplitudes of the magnetization 
on the coexistence line and the critical line with the result 

B = BT? = 1.245(7) ; d c = 0.978(2) , D c = 1.11(1) ; (32) 

where B = 0.945(5). The value for J c = T^ 1 = 0.454165 was taken from Ref. |T5|. 



3 Simulations at H = 



All our simulations were done on three-dimensional lattices with periodic boundary 



conditions. As in Ref. |T3| we have used the Wolff single cluster algorithm. The 



main part of the H = data was taken from lattices with linear extensions L = 
24, 36, 48, 72, 96 and 120. Between the measurements we performed 300-800 cluster 
updates to reduce the integrated autocorrelation time T int . Apart from the largest 
lattice (L = 120) where we made runs only at six couplings, we have generally 
scanned the neighbourhood of J c by runs at more than 30 points on each lattice, with 
special emphasis on the region 0.45414 < J < 0.45419. This enabled a comfortable 
reweighting analysis of the data. More details of these simulations are presented in 
Table [I]. 
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L 


J-range 


Nj 


iV meas [1000] 


Tint(t < 0) 


T in t(t ~ 0) 


T in t{t > 0) 


24 


0.440-0.4675 


35 


> 100 


1-3 


1-3 


1-3 


36 


0.440-0.4650 


43 


> 100 


1-4 


2-3 


2-10 


48 


0.442-0.4650 


55 


> 100 


1-5 


2-5 


4-13 


72 


0.4465-0.460 


41 


80-100 


1-4 


4-8 


7-21 


96 


0.450-0.4567 


33 


60-80 


2-10 


6-7 


7-35 


120 


0.452-0.4562 


6 


20 


2-4 


14 


12-23 



Table 1: Survey of the Monte Carlo simulations at if = for different lattices. Here 
Nj is the number of different couplings at which runs were performed; r int is the in- 
tegrated autocorrelation time for the energy and N meas the number of measurements 
per coupling in units of 1000. 

3.1 The Critical Point and the Binder Cumulant 

It is obvious that any determination of critical amplitudes relies crucially on the 
exact location of the critical point. Since we have produced a considerable amount 
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Figure 1: The Binder cumulant g r from Eq. @ as a function of the coupling J. 
The curves were obtained by reweighting the data. With increasing lattice size 
L = 24,36,48,72 and 96, the slope of the respective curve increases close to the 



critical point. The vertical dashed line denotes J c of Ref. ||16|| . 
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0.45415 0.45416 0.45417 0.45418 0.45419 0.45420 

Figure 2: The Binder cumulant g r in the close neighbourhood of the critical point. 
The figure is an enlargement of Fig. |l|. The dashed lines accompanying the solid 
lines show the jackknife error corridor. 



of data in the neighbourhood of the critical point it was natural to verify first the 
rather precise result of Ballesteros et al. |16| . We have done this by studying the 



Binder cumulant g r , which is directly a finite-size-scaling function 

g r = Q g (tL 1/u , IT ") . (33) 

The function Q g depends on the thermal scaling field and on possible irrelevant 
scaling fields. Here we have specified only the leading irrelevant scaling field pro- 
portional to L~ u , with uj > 0. At the critical point, t = 0, g r should therefore be 
independent of L apart from corrections due to these irrelevant scaling fields. In 
Fig. m we show our results for g r as obtained by reweighting the direct data. We 
observe, at least on the scale of Fig. [I], no deviation from the scaling hypothesis. 
However, after a blow-up of the close vicinity of the critical point, as shown in Fig. 
|2], we can see that the intersection points between curves from different lattices are 
not coinciding. The shift AJ of the crossing point from the infinite volume critical 
coupling J c can be estimated by expanding the scaling function Q g to lowest order 
in both variables. For two lattices with sizes L and V = bL one gets 

AJ L >" ex s(L,b) = l^jL— 1 '" ■ (34) 

In Fig. |3] we have plotted the J-values of the intersection points for each pair of 
lattices as a function of the variable s(L,b) of Eq. (|34]) . For uj we used the value 
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Figure 3: The coupling Ji P at the intersection point of g r (L) and g T (bL) for various 
combinations of L and h as a function of s(L, b), Eq. fl34"|). The filled (empty) symbols 
were calculated with v = 0.669 (0.673). The dashed (solid) lines are linear fits with 
(without) the L = 96 intersection points, denoted here by triangles. 



0.79(2) of Ref. []12| , and for v we have chosen the two values v = 0.669 and 0.673 as 
bounds of the probable iz-range. Of course, the intersection points are completely 
independent of v and u. Only the variable s(L, b) is changing when the exponents 
are changed. As can be seen in Fig. |3| also the extrapolation to the critical point J c 
for L — > (or s(L,b) — > 0) is unaffected by the choice of v. The same applies to 
a variation of u. Since the slope of g r (L = 96) close to the critical point is rather 
large, a small numerical uncertainty might shift the intersection points with the 
other curves considerably. We have therefore determined J c also by fits excluding 
the results from the largest lattice. Thus we arrive at the final estimate 

J c = 0.454167(4), (35) 



in full agreement with the result J c = 0.454165(4) of Ballesteros et al. [16|. In order 



to be consistent with our previous papers we use in the following again the value of 



Ref. 16 



In a similar manner one can determine from the same data the universal value 
g r {J c )- The difference of the g r - values at the intersection points to g r (Jc) is here 

Ag^ L ' ex So (L,b) = _ 1 . (36) 



8 



-1 .756 



-1 .758 



-1 . 760 



-1 . 762 - 



-1 . 764 



-1.766- 



. 



9 



P 



.05 



s (L,b) 

o ' 



.10 



Figure 4: The Binder cumulant g r ^ p at the intersection point for various combina- 
tions of L and b as a function of s Q (L,b), Eq. (|36|). The dashed (solid) lines are 
linear fits with (without) the L = 96 intersection points, denoted here by triangles. 



In Fig. f|we show the extrapolation of g r to the critical point value at s a (L, b) = 0. A 
variation of u in the range 0.77-0.81 leads only to a shift of 10~ 4 . The new variable 
s (L, b) is practically independent of u, the influence of v is not visible in Fig. £|. 
Comparing again extrapolations with and without the L = 96 points one obtains 



( ((J) 2 ) 2 ) 

g r (J c ) = -1.758(2) , or ^/(J e 

((f)- } 2 



1.242(2) 



(37) 



well in accord with the result of Ref . |TTJ (see also the long discussion in Ref . [ 17 ] ) . 



3.2 The Correlation Length 

In our H = simulations we have measured the correlation length using the second 
moment formula, Eq. (0). The finite-size-scaling equation for £ is 

£ = LQz(tL 1/u , L _£J ) , (38) 

and £/L = is a scaling function like g r , that is its value at the critical point is 
universal for L — > oo. In Fig. [5] we have plotted our correlation length data divided 
by L. Here formula ( |T4D has also been evaluated for J > J c or T < T c though in 
this region the data cannot be identified with the correlation length. We see again 
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Figure 5: The correlation length £ divided by L versus J for L = 24, 36, 48, 72 and 
96. The solid lines were calculated by reweighting the data. The dashed vertical 
line gives the postion of J c , the horizontal one the universal value, Eq. (|39f) . 



that all curves intersect at the previously determined critical point. A closer look 
into the neighbourhood of J c reveals however similar corrections to scaling as in the 
case of g r . The corresponding extrapolation of the variable s Q (L, b) to zero leads for 
i/L to 

i/L (J c ) = 0.593(2) . (39) 

This result confirms nicely the value £/L = 0.5927 from the preliminary simulations 
mentioned in Ref. [ |12|1 . 

Our data for the correlation length can also be used to find the critical amplitude 



£ + of Eq. (POD. To this end we use a method described in detail in Ref. |L 
We briefly repeat the main arguments assuming for simplicity that there are no 
corrections to scaling. An observable O with critical behaviour approaches for either 
positive or negative t and L — > oo the limiting form 

Ooo = a \t\- p , for |*1 — > , (40) 

where a is the critical amplitude and p the critical exponent. At finite L the 
observable satisfies a scaling relation 

0(t,L) = L p/u Q (x t ) , withx 4 = tL 1/v . (41) 

Here, Qo is the finite-size-scaling function of O. In order to ensure the correct 
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thermodynamic limit for fixed small \t\ we must have the relation 

Oo, = \t\- p lim \x t \ p Qo{x t ) • (42) 

The sign of xt is of course the same as that of t. It is clear then, that the function 

A (xt) = \x t \ p Qo{x t ) , (43) 

will converge asymptotically to the critical amplitude ao- Moreover, will be an 
extreme value of Ao(x t ). 

We have applied this method to the correlation length results. In Fig. |5| we show 
A^(xt) for the exponent v = 0.671 and various L-values. We notice that already 
at Xt ~ 4 a plateau is reached and essentially no corrections to scaling are visible. 
The marginal spread of the data in the plateau region leads only to a small error 
for the amplitude £ + . Since the scaling variable x t changes with v there is however 
a ^-dependence, which can also be expressed as a dependence on a, because of the 
hyperscaling relation 2 — a = dp. In fact, after evaluating for several ^-values, 
we find that £ + is rather exactly a linear function of the used a 

= 0.4957(20) + 0.67(12) a . (44) 

This can be seen in Fig. [7], where we compare the fit, Eq. to some directly 
determined £ + -values. 



. 5 





(x t ) 


1 










X§ x — m x-qx— * — e-OH 


8 

i 




□ 


96 


A 

- X 






72 


o 






48 






X 


B 








X 




A 


36 


o 




o 


24 


□ 














x t =tL 


X 




1 





5 10 

Figure 6: The amplitude function A%, Eq. ([43|), of the correlation length versus the 
scaling variable x 4 for v = 0.671 and L = 24, 36, 48, 72 and 96. The horizontal line 
indicates the £ + -value. 
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Figure 7: The critical amplitude £ + , Eq. Q2T|), of the correlation length versus —a. 
The data (circles) are determined from the amplitude function A^(x t ), the solid line 
is the linear fit (pj) . 

3.3 Specific Heat and Energy Density at T c 

As mentioned already in Section [2] both the energy density and the specific heat 
contain additional non-singular terms. This fact complicates of course the determi- 
nation of the critical amplitudes. We can however calculate the non-singular terms 
beforehand by a finite-size- scaling analysis directly at the critical point. For that 
purpose we have made further Monte Carlo runs at T c on 23 lattices with L = 8 to 
L = 160. In these runs we took between 500,000 and 200, 000 measurements each 
for L = 8 — 64 and on the larger lattices between 120,000 and 50,000. The data for 
the energy density and the specific heat are shown in Fig. |8] as a function of L up 
to L = 120. If one expands the scaling functions for e and C at T c in powers of 
one obtains 

e(L) = e ns + q 0e L^' v (l + q le L^ + . . .) , (45) 
C(L) = C ns + q 0C L a / v (l + q lo L-" + ...) . (46) 

We have fitted the first terms (up to q%) of these expansions to the data. In the 
case of the energy density we find no corrections to scaling, that is q\ e m 0, and 
only small corrections for the specific heat. Fits with different //-values cannot be 
distinguished in Fig. [8[ When we treat v as a free fit parameter we get v = 0.671(2). 
The quantity e ns exhibits no noticeable dependency on v or a and u. We find 

e ns = -0.98841(3) . (47) 
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Figure 8: The energy density (a) and the specific heat (b) versus L at the critical 
point. The dashed line shows e ns and the solid lines fits to Eqs. (|45| ) and ( [46l) for 
v = 0.671 and uo = 0.79. 



The situation is quite different in the case of the specific heat. Its non-singular 
part varies from about 50 for v = 0.669 to 16 at v — 0.675. The reason for this 
strong variation is that the exponent a = 2 — 2>v is close to zero, when v approaches 
2/3. Then the background term C ns develops a pole (~ 1/ct) which cancels a 
corresponding pole in the critical amplitude in such a way that the characteristic 
critical power behaviour (~ |£|~ a ) turns over into a logarithmic behaviour (~ In \t\). 
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Figure 9: The non-singular part C ns of the specific heat versus 1/a from fits to Eq. 
(f46|) (stars) with u = 0.79. The solid line is from Eq. fl54|). 

This mechanism for the emergence of the logarithmic singularity as a —>■ is well- 



known (see Refs. [15] and [19, P0|). We demonstrate it by assuming that 



C ns (a) 
A ± (a) 



c' 



i ~ns 
Wis T 

a 



(48) 

a± + afa + 0{a 2 ) . (49) 
If we insert these equations into Eq. (53) and expand |t| _a for small a we obtain 



C 



els + — + — + a T + O(a) (1 - a In |t| + . . .) 
a \ a 



_i_ "ns 



a 



+ af — Oq In |t| + 0(a) . 



(50) 
(51) 



Evidently the limit of C for a — > exists and has a logarithmic |t|-dependence, if 

a . (52) 



-do, and Oq 



the pole term vanishes, which requires (19| 

The ratio A + /A" is therefore close to 1 

A + /A~ = 1 + 0(a) . 



(53) 



In Fig. |9| we show the non-singular part C ns of the specific heat resulting from fits 
to Eq. (fHf) with = 0.79 and various values for a plotted versus 1/a. The \ 2 P er 
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degree of freedom in each fit is 0.83(1), preferring no particular a- value. We see 
that indeed C ns is linearly dependent on 1/q. A fit to the ansatz, Eq. (|48|), gives 



C n ,= 3.35(4) -™m, (54) 

a 

with an extremely small x 2 /Nf of the order of 10~ 4 . We conclude from this fact, 
that the pole term behaviour of C ns is not a numerical accident, but underlines the 
previous considerations. In order to study the influence of the correction exponent 
uj we have repeated the whole analysis of C(L) for the values uj = 0.77 and uj = 0.81, 
that is a standard deviation away from the central value 0.79. The x 2 /Nf for each 
single fit to Eq. flUp is again 0.83(1), the new values for C ns coincide within error 
bars with the values for u = 0.79, however the resultant linear fits in 1/a to Eq. 
(4£p at fixed uj, lead to slight changes (again with a \ 2 /Nf °f the order of 10~ 4 ) 

J 3.37(4) - 0.3165(5)/a for u = 0.77 . . 

ns ~ | 3.33(4) - 0.3184(5)/a for uj = 0.81 ' [ } 

mainly for the pole term parameter c£ s . 

In the following we shall use the results for C ns to analyze as well the specific heat 
data for T ^ T c . If not explicitly mentioned, the fit results have always been obtained 
for fixed uj = 0.79. We have repeated the following analysis also for uj = 0.77 and 
0.81 and shall comment on any noticeable changes due to uj. 



3.4 The Specific Heat and A + /A 



In Fig. [LOlwe have collected all our specific heat data at zero magnetic field for the L- 
values of Table [1]. We observe with increasing L a more and more pronounced peak 
close to J c . As already discussed in the introduction, we nevertheless expect a finite 
peak height even in the thermodynamic limit, since the singular part of C vanishes 
at the critical point for negative a. The peak (and not dip) behaviour implies also 
that the amplitude A ± /a must be negative, or that A ± is positive. The previous 
analysis of the non-singular contribution to C confirms this consideration: because 
c p ns is negative we have a positive value = a Q for the leading part of A ± . We 
have interpolated the data points by reweighting, apart from the L = 120 results. 
The respective curves are plotted in Fig. [11] as a function of t. Compared to Fig. 
|T0| we have therefore an exchange of the high (t > 0, J < 0) and low temperature 
(t < 0, J > 0) parts in the figures. In order to find the amplitudes A ± we have made 
the following ansatz including correction-to-scaling terms 

C = C ns + — \t\~ a (I + c±\tr + cft) . (56) 
a 

For a fit to the form (EB) the curves from the largest lattices were used in those 



t-ranges, which appear hatched in Fig. [11], that is for —0.0233 < t < —0.0045 and 
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Figure 10: The specific heat data for different L versus the coupling J. The dashed 
line indicates the position of the critical point. 



0.0048 < t < 0.0268. The non-singular part from Eq. (|54|) was then taken as an 
input to the fit, whereas the L = 120 data points served only as a check of the 
fit result. As an example we show in Fig. |11| the fit for a = —0.013. Fits with 
other small, negative a-values work as well and have the same x 2 P er degree of 
freedom, namely 1.03 . In Table |2| we present details of the fits for several a-values. 
The two correction-to-scaling contributions are always opposite in sign and cancel 
therefore to some extent, especially in the high temperature region. The amplitudes 
A ± are still a-dependent, though in our notation we have taken the anticipated pole 
behaviour already into account. We find that A + and A~ are nearly linear functions 
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Figure 11: The specific heat versus the reduced temperature t for L = 36, 48, 72, 96 
and 120 (stars). The solid lines were calculated by reweighting the data, the peak 
height increases with L. The line of long dashes is the fit from the ansatz, Eq. fl56|) , 
for a = —0.013 and u = 0.79. The hatched areas show the fit regions. 



a 


A+ 


4 


4 


A~ 




C 2 


-0.007 


0.3416(4) 


0.020(1) 


-0.041(1) 


0.3317(4) 


0.048(1) 


0.086(1) 


-0.013 


0.3636(6) 


0.022(1) 


-0.049(2) 


0.3445(6) 


0.085(1) 


0.161(2) 


-0.017 


0.3790(8) 


0.015(1) 


-0.041(3) 


0.3533(8) 


0.109(2) 


0.211(4) 


-0.019 


0.3870(9) 


0.010(2) 


-0.033(4) 


0.3578(9) 


0.120(2) 


0.237(5) 


-0.025 


0.4117(13) 


-0.016(3) 


0.006(6) 


0.3718(13) 


0.151(4) 


0.312(9) 



Table 2: The parameters of the fits to Eq. ( jSq) for uj = 0.79 and some selected 
a-values. The errors were obtained by Monte Carlo variation of the parameters of 
C ns in Eq. (|5^ 



of a. The a-dependence of the fit results for the amplitudes is shown in Fig. [12 
A parametrization of the amplitudes as suggested by Eqs. (|49|) and (|52| 



A ± = ao + afa + afa 2 



(57) 



works extremely well, as can be seen in Fig. 12, and confirms explicitly the cancella- 



tion of the pole terms as predicted in Eq. (|52f ) . If A + and A are independently fitted, 
that is with perhaps different ao, we get = 0.3176(12) and Oq = 0.3175(12). The 
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Figure 12: The amplitudes A + and A versus —a (squares). The filled circle is the 
value expected from C ns ; the lines are the parametrizations ( p8| ) and (|59|). 

final result is found by using Eq. fl57"D with fixed a = 0.3175 (the error in a = —c^ s 



is already included in the errors of the A ± - values, which are now parametrized) . We 
obtain 



A A 

A' 



= a - 3.308(36)a + 18.4(2.2)a 2 , (58) 
= a - 1.975(36)a + 7.8(2.2)a 2 . (59) 

At this point it is appropriate to discuss the influence of an ^-variation on A + and 
A~ . From Eq. (^) we know that a shift in u of size Au = 0.02 shifts the pole 
term parameter c? ns by about 0.3% and therefore we expect a shift of a by the same 
amount. In fact that is exactly what happens and it is the only effect, because 
the new parameters af and af coincide inside error bars with the values found for 
<jj = 0.79. All in all that results in a common shift of the A + and v4~-curves in Fig. 
|l~2l by again 0.3%. As a consequence the universal amplitude ratio A + / A~~ becomes 
essentially independent of uj. 



The universal ratio A + j A is sometimes given in terms of a function V(a) 

A + /A- = 1-Va . (60) 
Expanding the ratio in powers of a we arrive at the following relation for V(a) 



V 



-(1-4- 

a V A' 



Gh — Gh 



a 



+ 



a — an 



tt 



a 1 a 1 
a 



a 



a + 



(61) 



that is, V goes to a finite limit when a — > [P^, 52]. In fact, there is a phenomeno 

A+/A- 



logical relation M. 23 



1 -4a 



(62) 
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Figure 13: The universal ratio A + /A~ versus —a. The solid line is obtained from 
Eqs. fl58|) and (f)9|), the diamonds by direct calculation from Table |[ The other 
symbols denote results from the shuttle experiment (square) |7| ||, from Campostrini 
et al. (circles) [|1(], [LI], from Larin et al. (star) |24| and Kleinert et al. (plus) [ f25| . 



predicting V = 4. Evaluating Eqs. (|58|) and ( |59"D leads to 

= l-4.20(5)a + ... , 



(63) 



rather close to the relation (|62l). In Fig. [13] we show the ratio and compare it to 
former results from the shuttle experiment |7], B| as well as some analytical determi- 
nations [10, |TT| and [^4, 25| . We note that our ratio result is in complete accordance 
with all of the other ratio results. Obviously, they differ among each other simply 
and solely by assuming different a-values. This conclusion was already reached by 
Campostrini et al. flOfl, we can however directly confirm it with Eqs. and 



4 A + 1 A from the Equation of State 



The magnetic equation of state describes the critical behaviour of the magnetization 
in the vicinity of T c . As noted by Widom |19[ and Griffiths |22| already long ago 
the equation of state may be integrated to yield the scaling function for the free 
energy. From subsequent derivatives with respect to the temperature one obtaines 
then the specific heat and in particular an equation for the universal ratio A + /A~. 
Before we come to this relation we have to briefly discuss the equation of state. The 
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Widom-Griffiths form of the equation of state is given by 



y = fix) , (64) 

where 

y = h/M s , x = t/M l/f3 . (65) 
The variables i and h are the normalized reduced temperature and magnetic field 

t=(T-T e )/T , h = H/H , (66) 

associated with the usual normalization conditions 

/(0) = 1, and /(-1) = 0. (67) 

The reduced temperature i differs from t by a constant factor (t = [T c /To]t), because 
of the second condition in fl67]). The normalization constants can be expressed in 
terms of the critical amplitudes from Eq. (^) 

To = B-WT C = 1.18(2) , H = D c = 1.11(1) . (68) 

The numbers in the last equation have been obtained in Ref. [O] by assuming a 



special set |12] of critical exponents 

P = 0.3490(6) , v = 0.6723(11) , (69) 

which implies a ~ —0.017. The same is true for the equation of state, which was 
determined numerically in [13] from simulations with a non-zero magnetic field. 
Using this equation of state will therefore give A + / A~ for only that particular value 
of a. Varying a in the range [—0.0136,-0.0202], as suggested by the error of u, 
would result in a large variation of A + j A~ to begin with (see Fig. [13]). Insofar we 
consider the following calculation mainly as a test of the method. 

The results for the equation of state were parametrized in [|l^] by a combination 
of a small-x (low temperature) and a large-x (high temperature) ansatz. The small-x 
form x s (y) was inspired by perturbation theory p6| and incorporates the divergence 
of the susceptibility on the coexistence line (x = — 1; y = 0) due to the massless 
Goldstone modes 

x s (y) + 1 = (ci + d 3 )y + c 2 y 1/2 + d 2 y 3/2 . (70) 



The large-x form xi(y) was derived from Griffiths's analyticity condition [22 



xi(y) = ay 1 ^ + by^~ 2 ^ . (71) 

The parameter values are 

c x + d 3 = 0.352(30) , c 2 = 0.592(10) , (72) 

a = 1.2595(30) , b = -1.163(20) . (73) 
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Because of the normalization y(0) = 1 we have d 2 = 1 — {c\ + d^-\-C2)- The complete 
equation of state is obtained by interpolation of the low and high temperature parts 

x(y) = x s (y) P V ? - + xi(y) p V , (74) 
Vo + V P Vo + V p 

with p = 6 and y = 3.5. 

For negative a the universal ratio A + /A~ can be calculated from f(x) using the 



following formula |2y| 



-ffdxx"-* [f'(0)-f'(x) + f"(0)x] 

A' f(0)/(l - a) + /"(0)/a + $\ dx (-x)^ [/'(0) - /'(*) + /"(0)rr] ' 

The main contribution to both the nominator and the denominator is f"(0)/a. A 
more appropriate representation of A + / A~ is therefore 

A + l + [a/f"(0)]F N (?6) 



A- l + [a/f"(Q)]F D ' 



where 



f(0) 



1 — a 



F N = -f^-- / dxx a - 2 [f(0)-f(a;) + r(0)a:]+ / dxx a - 2 f(x), (77) 



o Jl 



F D = £22-+ / it(-ir J [f(0)-/'(i) + AO)i]. (78) 

1 - a y_i 

Let us denote the integrals in Eq. ( [77D by 7i and -Z2, the one in Eq. (|78|) by J3. 
To a good approximation we can calculate the integrals I\ and I3 as well as the 
derivatives from the low temperature equation (|7D|). In order to obtain I 2 we first 
rewrite the integral as 

r°° dr 

h = -/(l) + (2-a)/ dyy-fx a -\ (79) 

and evaluate the remaining integral from the interpolation formula (|74"|), using for 
/(l) the low temperature value 2.4448. For the derivatives we find 



/'(0) = 2^3-C!-d3-2c 2 ) 1 = 1.366 ±0.034, (80) 

/"(0) = [f(0)] 3 ((3/4)(c 1 + 4-l) + c 2 ) = 0.270 ±0.064, (81) 

and for the integrals 

h = 0.203 ± 0.02, I 2 = 1.749 ± 0.03, I 3 = 0.512 ± 0.02. (82) 
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The errors in the integrals were obtained by Monte Carlo variation of the initial pa- 
rameters in Eqs. ([72]) and (|73|). When this procedure is also applied to the complete 



expression (|76|) one obtains 

A + /A~ = 1.12 ±0.05 . (83) 

The first conclusion to be drawn from this result is that this method is not well suited 
for the calculation of the ratio, at least with the parametrization of the equation of 
state of Ref . |L3| . Though the result (|83|) is compatible with our directly determined 
ratio A + /A~(a = —0.017) = 1.073(3), the error is rather large. The main source 
of the error is evidently the inaccurate value of f"(0). That this quantity plays an 
important role is of course not unexpected, because A + and A~ are the amplitudes 
of the specific heat, which is again the second derivative of the free energy density. 
Our parametrization was not devised for that purpose, but for a correct description 
of the Goldstone effect near to x = — 1 and the limiting behaviour for x — > oo. That 
is why it led to a precise determination of R x and the constant Cf 

R x = lim x"</f(x) = 1.356(4) , c f = lim (1 +x)~ 2 f(x) = 2.85(7) . (84) 

x— >oo X— > — 1 

Campostrini et al. have used a different representation of the equation of state 



Pq , 11], based on Josephson's parametrization [29] of M,t and H in terms of the 
variables R and 9 and parametric functions. In order to fix these functions ap- 
proximately the authors utilized the results of an analysis of the high-temperature 
expansion of an improved lattice Hamiltonian. The values obtained for A + / A~ 
compare well with our direct determination and were already shown in Fig. [TB|. The 
corresponding equation of state differs however somewhat in the low and medium 
temperature regions from the data points from our non-zero field simulations [I3|. 



The question arises then whether the same data may be described as well in the 
schemes introduced by Campostrini et al. . Such alternative fits of the data have 
been carried out by two of us ][!(J. The \ 2 P er degree of freedom of these fits is 
generally high, in particular for scheme A of Ref. ||28|| . The fits according to scheme 
B are considerably better and lead to a ratio A + j A~ = 1.070(13), again compatible 
with our direct determination. The simultaneously calculated ratio Rc is however 
much larger (0.165-0.185) than expected from analytical calculations (0.123-0.130) 



31, 25]. We therefore do not pursue this method of calculation here in more detail. 



5 Simulations with H > 



We have performed additional simulations with a positive magnetic field H on the 
critical line to find the remaining critical amplitudes for the specific heat and the 
longitudinal and transverse correlation lengths. The linear extensions of the lattices 
we used were L = 36, 48, 72 and 96. These measurements were combined with those 
from Ref. Jl3| to cover the if-range appropriately. Some of the new data have 
already been used in Ref. [ 32| . In Table § we give more details of these simulations. 
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L 


ii-range 


1 y cu 


N meas [1000] 


N H 


N tot 


36 


0.0007-0.05 


50-100 


30-40 


25 


36 


48 


0.0001-0.03 


50-100 


30-40 


30 


39 


72 


0.0001-0.005 


60-300 


20 


15 


23 


96 


0.0001-0.0015 


60-80 


12-20 


8 


16 



Table 3: Survey of the new Monte Carlo simulations at T c on different lattices. N cu 
is the number of cluster updates between the measurements, N meas the number of 
measurements per if- value in units of 1000 and Nh the number of if- values at which 
new runs were performed. N to t is the total number of ii-values where we have data. 

5.1 The Specific Heat on the Critical Line 



In Fig. [14] we show our specific heat data as a function of the magnetic field H. 
Since there are no noticeable systematic finite size effects we can use these data to 
fit them to the ansatz 



C = C n 



A, 



H- a °(l + c h H u 



Here, C ns is the same non-singular term, which we have already determined in 
3.0 
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Figure 14: The specific heat at T c for L = 36, 48, 72 and 96 as a function of H. The 
line is the fit (B3) for a c = -0.0078 (u = 0.671) and u = 0.79. 
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a c 


a 


A c 


Ch 


X 2 /N f 


-0.00422 


-0.007 


0.2006(2) 


0.0203(1) 


1.09 


-0.00781 


-0.013 


0.2080(3) 


0.0344(2) 


1.09 


-0.01019 


-0.017 


0.2131(5) 


0.0423(4) 


1.10 


-0.01138 


-0.019 


0.2156(5) 


0.0458(4) 


1.10 


-0.01492 


-0.025 


0.2235(7) 


0.0546(8) 


1.11 



Table 4: The parameters of the fits to Eq. (|8~5"D for some selected a c -values at fixed 
f3 = 0.349 and u = 0.79. The errors were obtained by Monte Carlo variation of the 
parameters of C ns in Eq. (|88|) . 



Section [T3]as a function of a (or v) with the result (|54|). Because of the dependence 
of C on a c and v c the amplitudes A c and Ch depend on two critical exponents. The 
second exponent will however not introduce a sizeable variation in the amplitudes. 
We therefore treat the exponent (3 as fixed to the value (3 = 0.349 , in accord with 
our previous calculations. With the relations 

(35 = 2-(3-a, a c = -^ , a = a f~ P) , (86) 

2 — p — a l + a c 

the linear dependence of C ns on 1/a can be rewritten as one on l/a c 

°~ = c °- + 2^( 1 + i) (87) 

= 3.16(4)-°™. (88) 



We took this form of C ns as an input to the fits of C with Eq. (185]) . The if -range for 
the fits was 0.0001 < H < 0.05. We have convinced ourselves that smaller if-ranges 
(up to 0.02 or 0.03) lead inside the error bars to the same results for the amplitudes. 
In Table |] we present details of the fits for several a c -values, in Fig. ^ we show 
the amplitude A c as a function of a c . As in the case of the amplitudes A ± the pole 
of C ns in Eq. (|§S|) is compensated by the corresponding pole term in A c /a c . We 
have therefore parametrized the a c -dependence of A c in analogy to Eq. ( |5^) with 
the fixed value A c (a c = 0) = 0.1923 and find 

A c = 0.1923- 1.919(42)a c + 11.6(4.1)^ . (89) 

From Fig. |15| we see that this parametrization describes the data very well. Like in 
the study of the u- dependence of A^ in Section |3.4| we found changes of similar size 



for the amplitude A c due to a variation of u. They lead to an additional error of A c 
of size 0.0006 at a r = -0.00422, which decreases to 0.0004 at a r = -0.01492. 
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Figure 15: The amplitude A c versus — a c (squares) for uo = 0.79. The filled circle 
shows the value expected from C ns , the line is the parametrization fl59|) . 

5.2 The Correlation Lengths on the Critical Line 

The simulation results for the transverse and longitudinal correlation lengths are 
shown in Fig. [16] a) and b). For the transverse correlation length £r one can hardly 
detect finite size effects, whereas the longitudinal correlation length shows more 
fluctuations and a systematic deviation to higher ^-values, when one decreases the 
magnetic field H. The smaller the lattice, the earlier this behaviour sets in. In order 
to determine the amplitudes we have fitted our results to the following form 



£,T,L = &rH-"°(l + C TlL H'»">) 



(90) 





a 




c T 


a 


c L 


0.40350 


-0.007 


0.6709(14) 


0.024(13) 


0.3427(15) 


-0.258(33) 


0.40325 


-0.013 


0.6724(14) 


0.019(14) 


0.3435(15) 


-0.263(33) 


0.40307 


-0.017 


0.6735(14) 


0.015(14) 


0.3441(15) 


-0.266(33) 


0.40299 


-0.019 


0.6740(14) 


0.013(14) 


0.3443(15) 


-0.268(32) 


0.40274 


-0.025 


0.6755(14) 


0.008(14) 


0.3451(15) 


-0.273(32) 



Table 5: The parameters of the fits to Eq. ( ]9Ti| ) for some selected v c - values and 
uj = 0.79. The y 2 /Nf- values varied for £ T between 0.89 and 0.86, for £ L it was 0.67. 
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Figure 16: The correlation lengths £t (a) and £l (b) at T c for L = 36,48, 72 and 96 
as a function of H . The lines are the fits (pOf ) for v c = 0.40325 and uj = 0.79 . 



In the transverse case we used the reweighted data for L = 72 in the if -interval 
[0.0005,0.0025], for L = 48 in [0.002,0.02] and for L = 36 in [0.015,0.03]. From 
Table [| we see that the correction term is essentially zero. Correspondingly, there 
is no cu-dependence and a fit with cy = works just as well (even with the same 
X 2 /Nf), and leads to a slight increase in the amplitude value, which is of the order 
of the error given in Table |5|. The dependence of the amplitude £y on v c or a is 
linear but the slope is very small. In order to determine the longitudinal amplitude 
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££ we have fitted the reweighted data for L = 72 in the if-interval [0.0005,0.00175] 
together with those for L = 48 in [0.00175,0.01]. Here, the correction term is not 
zero, but the variation due to to is still negligible. The v c - or a-dependence is the 
same as for the ratio of the two correlations lengths is a fixed number 

C t /Cl = 1-957(10) , (91) 



independent of the critical exponents. It is well-known (see Refs. [[y| and [p3j p4| ) 
that at zero field on the coexistence line t < the longitudinal correlation function 
Gl is for large distances \r\ connected to the transverse one by 

G L (f,t) « ~(iV- l)[G T (f,t)/M] 2 , (92) 

where in our case N = 2. The relation is expected to hold also for small non-zero 
fields H near the phase boundary in the regime of exponential decay implying a 
factor 2 between the correlation lengths. It is remarkable, that we find approximately 
such a value for the ratio at t = 0. A similar observation has been made for the 



3d 0(4) model K5 



5.3 The Stiffness Constant on the Coexistence Line 

The stiffness constant p s (T) is related to the helicity modulus T |36| by 

Ps = T/T , (93) 



which can be measured in Monte Carlo simulations. This was done e. g. in Refs. |3j 
and |3*^| . Here we follow a different strategy, which we applied already in Ref. jT3| to 



find the magnetization on the coexistence line. The L or volume dependence of M 
at fixed J and fixed small H is described by the e-expansion of chiral perturbation 
theory (CPT) in terms of two low energy constants. One is the Goldstone-boson- 
decay constant F, the other the magnetization E of the continuum theory for H = 
and V — ► oo. The square of the constant F is proportional to the helicity modulus. 
In our notation, which is different from the one in CPT (see the remark in the last 
paragraph of Ref. P9]| ) we have 

T = F 2 /J , implying p s = F 2 . (94) 

The formulae, which are needed for the fits to determine the constants, are summa- 



rized in Ref. [0 and were taken from Ref. [59]. In Table |6| we list the results for the 
Goldstone-boson-decay constant F at various J-values. We performed simulations 
at H = 0.0001 on lattices with linear extensions L = 8, 10, 12, 16, 20, 24, 30, 36, 40, 48 
and 56. By construction the e-expansion is only applicable in a range where m T L<l. 
This condition translates into the equation 

H^=<(t)\ 05) 
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J= 1/T 


F 


AF 






0.462 


0.1993 


0.0096 


8,10,12 


36,40 


0.465 


0.2275 


0.0060 


8,10,12 


40,48 


0.470 


0.2596 


0.0050 


8,10,12 


40,48 


0.480 


0.3091 


0.0018 


8,10,12 


48 


0.500 


0.3795 


0.0114 


8,10,12 


48,56 


0.525 


0.4379 


0.0040 


8,10,12 


48,56 


0.550 


0.4755 


0.0028 


8,10,12 


56 



Table 6: The Goldstone-boson-decay constant F at various J- values from fits on 
data from lattices with L in the range [L min , L max }. 



and excludes the use of too large L-values. For each J we fitted different sets of 
data from lattices between [L min , L max ] and averaged the obtained F- values. The 
errors on F include the variations of these results. If we compare our F-values to 
the corresponding ones of Ref. we find generally somewhat lower numbers. This 
may be due to the fact that in Ref. [133] data from single lattices instead of sets of 



data from different lattices were fitted. The transverse correlation length £t on the 
coexistence line is now derived from the inverse of the stiffness constant or F~ 2 . It 




-0.20 



-0 . 15 



.10 



-0 . 05 



Figure 17: The inverse of the stiffness constant p s 1 = £t on the coexistence line 
from chiral perturbation theory. The solid line is the fit (|96|) with v = 0.671 and 
u = 0.79. 
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is plotted in Fig. [T7|. Here, we have not as many and as accurate data as in Fig. [16 
a). In order to determine the amplitude £y we fit our data points up to J = 0.525 
to the ansatz 

& = ^(-tr (i+c T Hr) . (96) 

Table [7| contains the fit parameters for different v or a-values. We observe, as for 
a linear dependence of the amplitude £^ on a with a very small slope. A change 
of uj by 0.02 leads only to a shift in £^ of a tenth of the error in Table [7|. 



V 


a 




Crp 


X 2 /N f 


0.6690 


-0.007 


1.680(52) 


-0.55(10) 


0.08 


0.6710 


-0.013 


1.665(52) 


-0.54(11) 


0.08 


0.6723 


-0.017 


1.655(51) 


-0.53(11) 


0.08 


0.6730 


-0.019 


1.650(51) 


-0.53(11) 


0.08 


0.6750 


-0.025 


1.636(51) 


-0.52(11) 


0.07 



Table 7: The parameters of the fits to Eq. fl9"6"| ) for several //-values and u = 0.79. 



6 The Universal Amplitude Ratios 

After having determined all the amplitudes which appear in Eqs. (|27j) to (|30"D we 
can calculate the corresponding universal ratios. Since the ratio Uq = A + /A~ has 
already been discussed in great detail we start with the ratio of the correlation 
lengths for H = 0. From Eq. ([44]) and Table [7| we find 

Of = t/^r = 0-293(9) , (97) 

independent of the used a-value. The e-expansion of this ratio was derived by 
Hohenberg et al. to O(e) and extended by Bervillier |3(J to 0(e 2 ) resulting in 



U{ = 0.27 and 0.33, respectively. Okabe and Ideura |4T| corrected the expansion 



of Bervillier (not the numerical value) and computed the ratio in 1/iV-expansion to 
U% = 0.140. The e-expansion results are comparable in size to our value in (p7[), the 
1 /iV-expansion result, however, seems to be too small. 

The ratios connecting the specific heat and correlation length amplitudes are 
related by 

R+ = RjU l/d Ut , (98) 

and they depend on the used a, mainly because of the specific heat amplitudes. In 
Table |8] we have listed the ratios and Rj. From the a-expansions ( [HQ and ( ]58"D 
we find 

R^ = 0.3382(14) - 0.717(96)a + 0.87(1. 13)a 2 . (99) 
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V 


a 


R t 


RJ 


Rc 


Ra 


Ql 


0.6690 


-0.007 


0.3432(15) 


1.163(36) 


0.118(4) 


0.0515(17) 


0. 


334(21) 


0.6710 


-0.013 


0.3476(18) 


1.167(36) 


0.125(4) 


0.0534(18) 


0. 


349(21) 


0.6723 


-0.017 


0.3505(21) 


1.170(36) 


0.130(4) 


0.0547(18) 


0. 


360(21) 


0.6730 


-0.019 


0.3520(22) 


1.171(36) 


0.133(5) 


0.0554(19) 


0. 


365(21) 


0.6750 


-0.025 


0.3563(27) 


1.176(36) 


0.142(5) 


0.0574(19) 


0. 


381(22) 



Table 8: The universal ratios from Eqs. (p8~D , (|29|) and ( |30l) as a function of the used 
exponents v and a. 



For Rj one can derive a similar formula representing the values of Table [8] 

Rj = 1.1580 - 0.696a + 0.97a 2 ± 0.036 . (100) 

There exist several theoretical estimates of R^ which compare well with our result: 
0.355(3) [a = -0.0146] and 0.361(4) @2|, both from high-temperature expan- 
sions; 0.36 g0| from the e-expansion, and 0.3597(10) g| and 0.3606(20) H from 



3d field theory. Apart from the first result, we could not relate a definite a- value to 
the respective estimate. The ratio Rj was calculated from the e-expansion |23| , |4~0| 
with the result 1.0(2) [jnj, well in accord with our value. 

The remaining universal ratios R x , Rc, Ra and Ql are all dependent on the 
amplitude C + of the susceptibility and/or the amplitudes B and d c (D c ) of the 
magnetization. We mentioned already that we had determined R x , B and d c in Ref. 



T3[| , although for fixed v = 0.6723. In the following we proceed as in Section |5.1| , 
that is we keep (3 fixed to 0.349 and assume in addition that the ^-dependencies 
of R x , B and d c are negligible. In Table [5] we present the ratios Rc and Ql as 
calculated from 

R c = A+R^B- 1 - 6 , Ql = (&/Z + )' l/u R x (dc/B) s - 1 /(l + l/5) , (101) 

and Ra directly from the definition in Eq. (0), using our newly determined am- 
plitudes A + ,A c ,£t an d £ + . We could not find any previous results for Ra and Q\ 
in the literature, however, the ratio Rc has been calculated theoretically in sev- 
eral ways. From Table |8] we see that Rc is increasing with decreasing a, which is 
due to the factor A + . In comparing our values to the analytical results we quote 
therefore the used a-values. The ratio Rc calculated from 3d field theory in Ref. 
3J is 0.123(3) [a = -0.01285], in Ref. 0.12428 [a = -0.01056]; from the high- 



temperature expansion in Ref. [[I fj one finds 0.127(6) [a = —0.0146]. The results are 
in full agreement with our calculation, though that of Ref. is somewhat higher 
than the other ones. The old e-expansion result 0.103 of Aharony and Hohenberg 



45 seems to be too small. 
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7 Conclusions 



We have calculated the major universal amplitude ratios of the three-dimensional 
0(2) model from Monte Carlo simulations. To reach this goal a large amount of 
computer time had to be spent on the cluster of alpha-workstations of the depart- 
ment of physics at the University of Bielefeld. Most of the computer time went into 
the production of reliable specific heat data for the direct determination of A + /A~. 
Initially we had hoped to improve the accuracy of the exponent a (or v) from these 
data. As it turned out, however, the specific heat data could be fitted to a whole 
range of a- values with the same \ 2 /Nf, extending even to a = 0. This raises the 
question, whether the experimental shuttle data are really fixing the a-value to 
exactly -0.01056, the same value as in 3d field theory expansions [|TJ. The positive 
aspect of the indifference of the fits to the specific heat data to a-variations was that 
we could study the numerical changes induced by these variations in the universal 
ratio A + /A~ and the background term C ns . As a result we were able to confirm 
the conjectured pole (in I /a) behaviour of the amplitudes and the background term 
and the mutual cancellation of the pole contributions. The same pole behaviour 
was observed for the specific heat amplitude on the critical line. The functional 
dependence of A + /A~ on the used a- value is in complete accordance with all other 
ratio results and not far from the phenomenological relation A + /A~ = 1 — 4a. We 
have also determined A + /A~ from the numerical equation of state, but we think the 
method relies too much on the chosen parametrization. 

In order to find the amplitude of the transverse correlation length on the coexis- 
tence line we used chiral perturbation theory. This enabled us to calculate the less 
known ratios Rj and U^. The latter is independent of the used a, like the ratio 
Ct/Cl on the critical line, which is remarkably close to 2 - a prediction expected for 
T < T c from the correlation functions close to the phase boundary. Our results for 
and Rc are in full agreement with the best theoretical estimates; Ra and Q\ 
are new and remain untested for the moment. 
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